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Abstract 

Microplasmas operated at atmospheric pressure show a number of peculiar dynamic phenomena. 
One of these phenomena is self-pulsing, which is characterized by intrinsic pulsing behavior of a DC 
driven plasma discharge. This work focuses on the numerical simulation of self-pulsing in a micro 
thin cathode discharge operated in atmospheric pressure argon. By means of a hybrid plasma 
model we show self-pulsing of the discharge in the expected MHz frequency range and described 
its actual origin. 
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FIG. 1: Sketch of the ^TCD with circuitry. Width of cathode and dielectric is 100 /mi, width of 
anode is 1 mm, hole diameter is 200 /jm. Applied voltage and series resistor are 1000 V and 25 
kfi, respectively. 

Microplasmas operated at atmospheric pressure have gained increasing attention since 
Schoenbach and co-workers did their pioneering work in the mid of the 1990th [1]. Mi- 
croplasmas show strong nonequilibrium behavior although they are collisional [2] . Electrons 
are hot (up to a few eV) whereas ions and atoms (and/or molecules) of the neutral gas are 
cold (about room temperature). Microplasmas are characterized by high electron densities 
(up to a few 10 17 cm -3 ) and high power densities (up to some 100 kW/cm 3 ). Due to this 
extraordinary properties microplasmas have found widespread of technological and biomed- 
ical applications. Examples can be found in the area of surface treatment, sterilization, and 
lighting [3, 4]. 

Microplasmas show unexpected dynamical behavior. One particular phenomenon that is 
experimentally observed in DC driven micro discharges is self-pulsing, which occurs as quasi 
periodic peaks in the discharge current. It is argued that this phenomenon can be understood 
as the periodic transition from the Townsend regime to the abnormal glow and vice versa. 
For a micro hollow cathode discharge operated at gas pressures between 50 and 100 Torr 
Rousseau and co-workers experimentally found the pulsing frequency ranging from a few 
tens to a few hundreds of kilohertz [5, 6]. For the micro thin cathode discharge (/iTCD) 
proposed by Czarnetzki and co-workers the self-pulsing phenomenon is in the megahertz 
range [7]. 

Before the advent of microplasmas Petrovic and Phelps discussed certain oscillations of 
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low current electrical discharges between plane-parallel electrodes [8-10]. They propose 
lumped circuit models in order to describe the interaction of the nonlinear negative differen- 
tial resistance "behavior" of the discharge with the external circuit. To study the self-pulsing 
phenomenon of micro discharges Hsu and Graves did pioneering work [11]. They also pro- 
pose a lumped circuit model for the self-pulsing regime of micro hollow cathode discharges. 
Chabert et al. acted on the suggestions of Hsu and Graves and set up an improved lumped 
circuit model in order to carefully study the dynamics of the system [12]. Kolobov et al. 
described the transition from a Townsend discharge to a normal discharge using a two- 
dimensional numerical model [13]. Their numerical model is based on a fluid description 
of electron and ion transport coupled with Poisson's equation, with the ionization source 
depending on the local field strength or provided by a Monte Carlo simulation of the fast 
electrons. The model is applied to an argon discharge, for a product of pressure and gap 
length in the 1-10 Torr cm range (which holds also for microplasmas) . In 2003 Arslanbekov 
and Kolobov re-investigated the phenomenon by coupling two commercial simulation tools 
[14]. Low-current self-generated oscillations in a rectangular hollow cathode discharge in he- 
lium gas were investigated using a two-dimensional self-consistent hybrid model by Donko 
[15]. The model combines Monte Carlo simulation of the motion of fast electrons and a fluid 
description of slow electrons and positive ions. The low-frequency oscillations of approx. 
20 kilohertz were found to arise as an effect of the interaction of the gas discharge and the 
external electric circuit consisting of a stable voltage source, a series resistor and a capacitor 
formed by the discharge electrodes. 

In this work we study self-pulsing of the DC driven /zTCD. The cross section of the 
/xTCD is schematically depicted in Fig. 1. The difference between the micro hollow cathode 
discharge described by Schoenbach [1] and the /iTCD investigated here is that the cathode 
is significantly thinner than the anode. The cathode has a thickness of 100 /im whereas the 
anode is 1000 /im thick. The dielectric has a thickness of 100 /im. The hole has a diameter 
of 200 /im. 

In order to study the dynamics of the discharge we apply the two-dimensional hybrid 
plasma model nonPDPSIM by Kushner and co- works, which is described in detail in [17]. 
Here we briefly discuss the relevant physical equations. The electrical potential is calculated 
in electrostatic approximation giving Poisson's equation. Coupled to drift- diffusion equa- 
tions charged species and a surface balance equation self-consistency is reached. The set 
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FIG. 2: The discharge current (mA) over time (ns). The vertical lines at 170 ns, 216 ns, 223 ns, 
230 ns, 234 ns, 240 ns, 252 ns, and 300 ns correspond to electron density plots in Fig. 4 

of equations is simultaneously integrated in time using an implicit iteration scheme. Up- 
dates of densities and mean velocities of charged species are followed by an implicit update 
of electron temperature by solving the electron energy balance equation. To capture the 
non-M axwellian behavior of electrons electron transport coefficients and rate coefficients are 
obtained by solving a zero-dimensional Boltzmann equation for electrons. Secondary elec- 
trons are treated in frame of a Monte-Carlo simulation. The discharge is sustained in argon 
at atmospheric pressure. The species in the model are electrons, Ar(3s), Ar(4s), Ar(4p), 
Ar + , Arj^ and ArJ. The reaction mechanism is summarized in Ref. [17]. Due to its the 
small thickness the cathode is easily heated by the plasma leading to thermionic electron 
emission from the cathode surface. This significantly contributes to power deposition into 
the discharge. However, this effect can clearly be separated from the self-pulsing and is 
therefor not accounted for in this afford. 

The simulation addresses a time interval of 900 ns and is performed on a computational 
domain of 500 fim width and 1 mm height. A DC driving voltage of -1000 V is applied 
to the cathode via a series resistor of 25 kQ. Electron emission is described applying a 7 
coefficient of 0.15 for the cathode surface facing the plasma region. 
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FIG. 3: The cathode voltage (V) over time (ns). The vertical lines at 170 ns, 216 ns, 223 ns, 
230 ns, 234 ns, 240 ns, 252 ns, and 300 ns correspond to the electron density plots in Fig. 4 

The electron density as a function of time in the interval of 170 ns to 300 ns is presented 
in Fig. 4. At the instant time of 170 ns the simulation has reached a quasi steady state with 
a flat maximum of about 2.5 x 10 15 cm -3 at the edge of the cathode. A potential plot (not 
shown here) shows a equipotential area in the surrounding of this maximum, indicating a 
glow mode plasma. We like to call this area core in the following. In the following 60 ns a 
second maximum in the same order of magnitude is formed in the vicinity of the dielectric 
(Fig. 4b, c). In contrast to the core this maximum has a peak shape. It appears that this 
density peak is the nucleus of a pulse. In Fig. 4e the peak starts to advance towards the 
cathode and merges finally with the core. Another density peak forms approximately 40 ns 
later and the process is repeated. Interestingly the density peak launches from different 
positions along the dielectric resulting in a non-constant pulsing frequency. 

Fig. 2 illustrates the electron and ion density distribution at the instant of time of 
233 ns, when the simulation reached a quasi steady state. The discharge shows a glow-like 
character with a flat maximum of about 2.5 x 10 15 cm~ 3 at the edge of the cathode. We obtain 
another maximum of the electron density of the same order in the vicinity of the dielectric. 
In contrast to the flat maximum at the cathode edge this maximum shows a peak shape. 
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FIG. 4: Electron density in plasma region as function of time on log scale 



Looking at the dynamics of the plasma this maxima appears to be the beginning of a pulse. 
We focus on the phenomenon in Fig. 4. Here the temporal evolution of the inner maximum 
is given. The peak expands and moves in a beam towards the main maxima over a period of 
approx. 120 ns and eventually vanishes. Another peak is formed approx. 30 ns later and the 
whole process starts again. Interestingly in each period the "beam" launches from different 
positions along the dielectric and the pulsing frequency is not constant. Consequently, the 
discharge shows chaotic behavior in the given parameter regime. Fig. 3 shows the discharge 
current over a period of 1000 ns. Here the self-pulsing behavior can clearly be seen as well. 
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